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Abstract. Context. High contrast imaging at optical wavelengths is limited by the modest correction of conventional near-IR 
optimized AO systems. We take advantage of new fast and low-readout-noise detectors to explore the potential of fast imaging 
coupled to post-processing techniques to detect faint companions of stars at small angular separations. 

Aims. We have focused on /-band direct imaging of the previously detected brown dwarf binary HD 130948 BC, attempting to 
spatially resolve the L2+L2 system considered as a benchmark for the determination of substellar objects dynamical masses. 
Methods. We used the Lucky-Imaging instrument FastCam at the 2.5-m Nordic Telescope to obtain quasi diffraction-limited 
images of HD 130948 with ~0.1" resolution. In order to improve the detectability of the faint binary in the vicinity of a bright 
(7=5.19+0.03) solar-type star, we implemented a post-processing technique based on wavelet transform filtering of the image 
which allows us to strongly enhance the presence of point-like sources in regions where the primary halo generally dominates. 
Results. We detect for the first time the binary brown dwarf HD 130948 BC in the optical band / with a SNR-9 at 2.561"±0.007" 
(46.5 AU) from HD 130948 A and confirm in two independent dataset (May 29 and July 25 2008) that the object is real, 
as opposed to time-varying residual speckles. We do not resolve the binary, which can be explained by astrometric results 
posterior to our observations that predict a separation below the telescope resolution. We reach at this distance a contrast of 
A/=11.30±0.11, and estimate a combined magnitude for this binary to /=16.49±0.11 and a I-J colour 3.29±0.13. At 1", we 
reach a detectability 10.5 mag fainter than the primary after image post-processing. 

Conclusions: We obtain on-sky validation of a technique based on speckle imaging and wavelet-transform post-processing, 
which improves the high contrast capabilities of speckle imaging. The I-J colour measured for the BD companion is slightly 
bluer than, but still consistent with what typically found for L2 dwarfs (~3.4-3.6). 
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1. Introduction 

A direct determination of dynamical masses of very low 
mass (VLM) objects is essential to calibrate the mass- 
luminosity relationship. This is particularly relevant for un- 
derstanding brown dwarfs (BDs) evolution. Dynamical masses 
can be determined by observing close multipl e BD systems 



(Zap atero Osorio etall 2004c IBouv et al.l l2004t Istassun etal. 
2006I iDupuv et alJ l2009h . BD close binaries with orbital 



periods <10yr represent a precious sample for a model- 
independent mass determination within a realistic time base- 
line. Observationally, this requires to spatially resolve the bi- 
nary, which also permits us to obtain a direct measurement of 
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the flux of each component. Since BD systems are also de- 
tected as close companions to bright main-sequence stars, an- 
other difficulty resides in the str ong contra s t need ed to detect 
them (cf. the case HR7672B in [Liu et all J2002h ). on top of 
the detectability issue due to their intrinsic low luminosity. So 
far, the sample of such companion B P binaries is limited to 
a few number dBurgasser et al.ll2005l) . mostly characterized in 
the near-IR with the help of 8-10m class telescopes. 
Optical data are necessary for a full characterization of the 
spectral energy distribution, key to the determination of ef- 
fective temperatures and bolometric luminosity. In the vis- 
ible domain, clos e binaries can b e spatially resolved using 
speckle imaging (ILaw et al.ll2006al) . a technique that delivers 
diffraction-limited optical counterpart to AO-assisted infrared 
images. The question of high contrast in speckle imaging has 
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been investigated in the past bv iBoccaletti et aL (2001) using 
the "dark speckles" method as an additional stage of clean- 
ing to improve the detectability of faint companions. Coupled 
to the adaptive optics system ADONIS and a Lyot stellar 
coronograph, these authors obtained K-band contrasts of 1.5- 
4.5xl0~ 3 {Am K ~6-l) at 0.5-0.9". 

In this paper, we have focused on t he brown dwarf bi nary 
HD 130948 BC, originally reported bv lPotter et ail d2002h . As 
part of a restricted sample of BD binaries companion to a 
solar-type star, HD 130948 BC is a unique benchmark for the 
study of the mass, luminosity and age of L-type substellar 
objects. Lying at -45 AU from HD 130948 A, the determina- 
tion of the clo se BD binary orbi t based on infrared AO im- 
ages permitted iDupuv et al. (2009) to derive a total mass es- 
timated to 0.109+0.03 M Q . We intent to use speckle imaging 
techniques to conduct a first optical ground-based high resolu- 
tion study of this object in order to extend its physical charac- 
terization to shorter wavelengths, and at the same time using 
it as an observational testbench for diffraction-limited imaging 
from the ground at visible wavelengths. Here, we report first I- 
band high contrast speckle imaging of the HD 130948 system. 
The detection of the BD companion is improved using image 
post-processing. This emphasizes how small-medium size tele- 
scopes can still be exploited for high angular and high contrast 
imaging. In § 2 are presented our observations and the data re- 
duction procedure. In § 3 we present our imaging, astrometric 
and photometric results on HD 130948 BC, briefly discussing 
our derived I - J colour. 



2. Observations and Data Reduction 

2.1. Imaging 

HD 130948 was observed on the nights of 2008 May 29 and 
2008 July 25 using the FastCam instrument (lOscoz et al . 2008) 



installed at the 2.5-m Nordic Telescope (NOT) at the Roque de 
los Muchachos Observatory, Spain. Observing conditions dur- 
ing these nights were good with an average seeing of 0.5" in the 
/-band and clear weather. In brief, FastCam is an optical imager 
based on a conventional low-noise CCD camera from Andor 
Technology that allows us to record speckle-featuring unsat- 
urated images at a rate of several tens of frames per second. 
Each frame captures a different pattern introduced by the atmo- 
spheric turbulence, where each speckle represents a diffraction- 
limited image of the source of interest. 

FastCam was installed on the Cassegrain focus of the NOT. 
The pixel scale, determined from astrometric measurements in 
the Ml 5 globular cluster, is estimated at 31.17+0.03 mas/pixel. 
The field-of-view accessible with the 512x512 CCD array is 
16"xl6". The raw data acquired by FastCam are composed 
of cubes of 1000 images each, with an individual integration 
time of 10-50ms. A total of 100 cubes (i.e. 10 5 frames) and 
50 cubes was acquired for HD 130948 in the nights of July 25 
and May 29, respectively, as well as flat and dark images for 
standard image correction. After finding the brightest speckle 
in each image of the cube, we can select an arbitrary percent of 
images based on the brightest speckles sorted within the serie. 
Hence, are selected those frames with the highest concentration 



of energy in a given speckle. A subsequent shift-and-add pro- 
cess of the previously selected frames products an image where 
the bright diffraction-limited core is surrounded by a fainter 
seeing halo. This technique is also referred as Lucky Imaging 
dLaw et al.ll2006 b). The resulting angular resolution depends 
then on the percent of selected frames and, among others, on 
the natural seeing during the night and the integration time for 
each individual frame. The percent of selected frames results 
from a trade-off between a sufficiently high integration time 
and a good angular resolution. In our case, the individual in- 
tegration time and the percent of selected frames were fixed, 
respectively, to 30 ms and 30%, leading to an effective total 
integration time of 900 s. Thus, FastCam delivered ~0.1" res- 
olution images, close to the diffraction-limit 1 .22/l/D of 84 mas 
of the NOTat0.8//m. 

As in any imaging system aiming at compensating the effect 
of the atmospheric turbulence, the presence of the seeing halo 
in the images is a strong limitation to high contrast capabili- 
ties. Since we are not able, with the current observational tech- 
nique, to artificially reduce the starlight contribution as it can 
be done using a coronagraph, the flux contrast at a given dis- 
tance from the star is improved by applying an additional post- 
processing stage to the image, which helps to suppress most 
of the seeing halo, unveiling fainter sources in the immediate 
vicinity of the star. We initially tried the classical solution of 
PSF-reference subtraction. However the strong time-variability 
of Lucky-Imaging PSF profiles at optical wavelengths prevents 
us from applying an efficient subtraction that would improve 
the detectability. Thus, we implemented a different approach 
based on the post-processing of the shift-and-add images in 
order to enhance particular spatial frequencies in the image. 
For this purpose, we have tested different image filtering algo- 
rithms, which are: 

★ the subtraction of a median box filtering. 

★ the implementation of a standard unsharp mask (i.e. I=I-I*g, 
where / is the input image, "*" denotes convolution and g is a 
Gaussian kernel. 

★ the implementation of the wavelet transform of the shift-and- 
add image. All these algorithms, comparable to edges enhance- 
ment techniques, allow us to suppress a continuous offset (e.g. 
sky background, detector bias...) or an extended and diffused 
structure (e.g. a PSF halo) at a different spatial scale than the 
object of interest, hence favoring the detection of fainter point- 
like sources at higher spatial frequencies. 

The best results were obtained with the unsharp mask filtering 
and the wavelet decomposition. In their numerical implemen- 
tation, those two algorithms are very similar, but the wavelet 
algorithm gives the advantage of a multi-resolution approach. 
This last solution was finally selected and applied to our two 
different dataset of HD 130948, and is described in more de- 
tails in AppendixlAl 



2.2. Photometry 

We used an absolute measurement of the HD 130948 A pho- 
tometry to derive the / magnitude of the companion from our 
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images (see Sect. 13.41 1. In order to determine the / magnitude 
of HD 130948 A, we carried out a dedicated photometric mea- 
surement obtaining /-band images of HD 130948 on 6 April 
2010 using the CAMELOT instrument mounted on the IAC80 
telescope. This optical camera consists of a 2kx2k CCD de- 
tector with a 0.304"/pix plate scale providing a field of view 
of 10.4x10.4 arcmin. The same / filter can be physically inter- 
changed in FastCam and CAMELOT and match the Johnson- 
Cousins system. The spectral responses of the each of the two 
detectors were obtained from available technical data and the 
differential effect on a solar-type star spectrum was estimated 
to a residual magnitude of 2.9milli-magnitude. We observed a 
serie of 10 images with individual exposure time of 5s. Raw 
data were reduced using routines within the IRAF environ- 
ment. Bias-images were subtracted using the overscan region 
and zero exposure time images, and flat-field correction was 
applied using dome flats. We performed aperture photometry 
using routines from the DAOPHOT package. We adopted an 
aperture of 5 FWHM (depending on the seeing conditions, this 
number varies bet ween 6 . 4" to 8 .7") for our photometric stan- 
dard stars from lLandoltl dl992l) . This is a standard aperture 
recommended by several photometric manuals to avoid con- 
tamination. Since the HD 130948 images were defocused to 
avoid saturation, we adopted a larger aperture (10 FWHM) and 
hence, in order to correct for the difference, we estimated the 
aperture correction from 5 to 10 FWHM using bright and iso- 
lated Landolt standard stars. We obtained an aperture correc- 
tion of 0.033 +/- 0.005 mag, which was included in the photom- 
etry of our target. Weather conditions during our observations 
were photometric as assessed by observing photometric stan- 
dard all through the night, while average seeing ranged from 
1.3 to 2". In order to transform our instrumental magnitudes 
into apparent magnitudes, we observed four different Landolt 
standard star fields (each of them containing 3-6 standard stars) 
and repeated them along the night. We obtained 14 different im- 
ages at 7 different pointings covering a range of airmasses from 
1.1 to 2. 1 . We perform a linear fit to our data to obtain the zero 
points and the extinction coefficient, following the equation 
i-I-aQ+kxairmass, where i and / are, respectively, the instru- 
mental magnitude and the apparent magnitude of the Landolt 
stars, ao is the zero point, k is the extinction coefficient. We 
obtained a =2.68+0.019 (zmag=25) and £=0.17+0.013. The 
error bars in the calibration were obtained from the estimated 
errors of the coefficients in the linear fit. Eventually, our correc- 
tion from instrumental to apparent magnitudes for HD 130948 
was finally i - 1= 2.860+0.032 (-0.033+0.005 aperture correc- 
tion ) = 2.827+0.032. The final error bar includes both the error 
in the calibration and in the instrumental magnitude. 

3. Results 

3.1. Detection of the BD companion 

The panels of Fig.[T]show the imaging results from our obser- 
vations . In the upper-left c orner is displayed the original detec- 
tion bv lPotter et alJ (I2002I) in which the brown dwarf binary is 



HD 130948 obtained with the data of July 25. In all the images, 
North is up, East is left. The FastCam images have been rotated 
by 90° with respect to the original position of the detector on 
sky. The average full-width-at-half-maximum (FWHM) is 131 
mas, with a slight elongation in the East direction (148 mas) 
against the North direction (114 mas). We attribute this effect 
to atmospheric dispersion since no ADC (atmospheric disper- 
sion compensator) is available in the current FastCam+NOT 
configuration. The white arrow indicates the expected position 
of the BD companion with respect to HD 130948 A, while the 
inset shows, on a different intensity scale, the core of the PSF 
The bottom part of Fig. Q] shows the HD 130948 system ob- 
served in July (left) and May 2008 (right) after the image fil- 
tering step. The brown dwarf is detected with a signal-to-noise 
ratio (SNR) ~9 and an average FWHM~1 10 mas, at a position 
very consistent with earlier image (see next section). The BD 
binary HD 1 30948 B - HD 1 30948 C is unresolv ed in our im- 
age. A ccording to the orbital solution derived by Dupuv et al 



resolved in the //-band with the 8-m Gemini-North telescope. 
The upper-right corner shows the direct shift-and-add image of 



(2009), the BD separation at epoch 2008.56 is estimated to 
30 mas, below the theoretical diffraction-limit resolution of the 
NOT. The physical separation at the same epoch is 0.63 AU 
assuming c/=18.17pc given by the same authors. Because the 
filtering process leaves speckles in the image that may mimic 
the presence of a companion, the detection is doubled-checked 
using a procedure that splits the 100 cubes of the July dataset 
into three separate datasets, which undergo the same post- 
processing step. In this way, speckle features simply resulting 
from the image processing are unlike to be repeated in the dif- 
ferent sub-images, while a real companion remains detectable, 
although with a lower SNR. In our particular case, we benefit 
also from the additional dataset of May 29 (50 cubes) in which 
HD 130948 BC is searched for. In both cases, the BD compan- 
ion binary is detected at the same location with respect to the 
G2V star. In the following part of this paper, the shift-and-add 
image of the best 30% selected frames and the filtered image 
will be referred, respectively, as Image 1 and Image 2. 



3.2. Detectability and contrast curve 

In the Fig. [2] we present the the 3-cr detectability curve as 
a function of the distance to the central star for Image 1 and 
Image 2. These curves were obtained by estimating, in the cor- 
responding images, the radial standard deviation profile in a 
9x9 pixels slipping box averaged in azimuth over 22 angular 
positions in clean areas of the image (i.e. not affected by the 
vertical and horizontal spikes visible in Fig. [TJ and separated 
by 5°. The gain is of the order of 2 mag when implementing the 
post-processing step. This result, interpreted in the perspective 
of the images of Fig. [T] suggests that the presence of the bright 
PSF halo in the direct image significantly contributes to the 
background noise at various spatial frequencies, degrading the 
detectability of a faint companion. On the contrary, in the fil- 
tered image the background noise associated to the bright halo 
is mostly filtered out around the spatial frequency constrained 
by the wavelet kernel. In other words, the wavelet filtering step 
strongly improves the signal-to-noise ratio of the object detec- 
tion in absence of any external hardware system aiming at re- 
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Fig. 1. Upper-Left: Original //-band detection of HD 130948 BC in Potter et alj (12002b . - Upper-Right: Direct lucky imaging 
image resulting from the best 30% frames over a 10 5 serie. The white arrow indicates the anticipated position of the companion. 
The inset shows, on a different intensity scale, the core of the PSF. - Bottom-Left: Post-processed image from July 2008 revealing 
the BD binary companion, unresolved with the NOT at 0.8 pm. - Bottom-Right: Post-processed image from May 2008 obtained 
with 5xl0 4 frames. All four images have the same size. North is up, East is left. The effective total integration time to obtain this 
image is 900 s for the July image. 



ducing the PSF halo (e.g. with a coronograph). In Fig. [2] was 
added a point showing the detection of HD 130948 BC (see 
next section for contrast estimation), which helps to better per- 
ceive the gain set by the post-processing step. 

3.3. Astrometry 

To determine the separation between the primary and the BD 
binary, we first used the photo-center of Image 1 (or barycenter 
in a 15x15 pixels aperture centered on the primary) as the astro- 
metric reference since we benefit from an excellent signal-to- 
noise in the unsaturated image of the primary. The position of 
the photo-center coincides with the brightest pixel within 0.09 
pixels. The same result was obtained using a two-dimensional 
Gaussian fit of the primary PSF, which is explained by the 
high SNR obtained on this bright object. At the contrary, we 
only implemented the Gaussian fit technique to determine the 
position of the faint companion in the filtered image, which 
proved to be more accurate than the photo-center technique 
for lower SNR. We first checked, by simulating several ar- 
tificial companions at similar distance and with comparable 



brightness to HD 1 30948 BC, that the presence of such a faint 
companion has negligible impact on the position of the pri- 
mary photo-center. We also verified that the position of the 
faint companion in the filtered image (Image 2) is possibly al- 
tered by the process of spatial filtering by less than 0.2 pix- 
els (or 0.007") in average. The advantage of measuring the 
BD position in the filtered image is to benefit from an en- 
hanced detection, removing the limitation of a poor signal- 
to-noise ratio. The two-dimensional Gaussian fit provided an 
estimation of the relative positions of the two components 
of 2.561"+0.007" at epoch 2008.56 (46.5 AU at 18.17pc), 
slightly closer compared to the separa ti on of 2.64"+0.01" 
at epoch 2001.15 found by [Potter et all d2002l) . The relative 
orientation of HD 130948 BC with respect to the primary is 
PA. = 102.7°+0.15°, to be compared with PA. = 104.5°+0.5° 
foun d by the previous a uthors. Assuming the distance revised 



by iDupuv et al 



J2009I) to d=18.17+0.11pc, these measure- 
ments may indicate an orbital motion around the primary, with 
a projected physical separation decreased by 1.43+0.22 AU in 
7.4 yr and a PA. change by 1.8°±0.55° in the same period. 
These values are also in agreement with the astrometric data 
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Arcsec 

Fig. 2. 3-cr detectability curve estimated in the two cases of 
Image 1 (dashed line) and Image 2 (continuous line). These 
data were obtained with 100 cubes and a 30% best-frames 
selection per cube. The filled triangle shows the position of 
HD 1 30948 C (see text for details). 



provided by Geifile r et al. I d2008h in their mid-infrared observa- 
tion of the HD 130948 system at epoch 2006.51. This change 
corresponds to an approximate orbital velocity projected onto 
the sky of ~1 km.s -1 , to be compared with an unprojected or- 
bital velocity of ~4.6km.s~ 1 assuming a circular orbit of ra- 
dius 46.5 AU and a stellar mass of 1.11 M for HD 130948 A. 
The position of the object has been also measured in the image 
from May observations. Although the estimation is less pre- 
cise (+0.015" and ±0.3°), we find consistent values between 
the two epochs with 2.56" and 102.5° at epoch 2008.41. 



us to directly compare, in the processed image, the peak values 
between the synthesized companion and the original detection. 
The significance of this approach is justified by the small width 
of the field-of-view (±3"), which is within the typical isopla- 
natic patch. As a consequence, the speckle pattern in the raw 
short-exposure is not expected to show significant spatial vari- 
ations within +3", with comparable implications on the PSF 
of Image 1 . Hence, for each contrast value in the image model, 
we benefit from 25 independent experimental contrasts from 
which are extracted mean and rms values. The statistical dis- 
persion express possible PSF distortion effects in the field-of- 
view affecting the peak value of the companions (see Fig. [3). 
From this empirical relationship, the contrast of HD 1 30948 BC 
is estimated by interpolation to the measured peak value. As a 
second approach we have also applied the reverse reasoning, 
which consists in tuning to an arbitrary precision the fake com- 
panion contrasts in the image model to match the experimen- 
tal peak value of HD 13948 BC in the post-processed image. 
We then derive a new estimate of the mean contrast and corre- 
sponding dispersion. Eventually, we apply the same procedures 
using as input measurements the fluxes summed in the aperture 
described above. 

We are provided in the end with four different estimates 
of HD 1 30948 BC contrast with a conservative value of 
3 .0+0.3 x 1 0~ 5 ( A7= 1 1 . 30+0. 1 1 ) with respect to the 7=5 . 1 9 pri- 
mary. Note that the peak value method provides a better accu- 
racy (+0.21 xl0~ 5 ) than the aperture photometry (±0.32xl0~ 5 ), 
in part because in the last case we do not weight the con- 
tribution of each aperture pixel as this can be done with op- 
timal extraction techniques. The conversion into magnitudes 
results into an unresolved BD magnitude 7=16.49+0.11 (or 
7=17.24 per component assuming equal brightness for the 
two c omponents). The c ombined MKO photometry of the BD 
from iDupuv et alJ d2009h results in the 7=13.2+0.08 magni- 
tude, leading to a colour index 7-7=3.29+0.13 for the unre- 
solved binary. We have also verified that the possible elonga- 
tion by ~30 mas of the unresolved binary has a negligible effect 
(<5 milli-mag) on the measured photometry. 



3.4. Photometry of HD 1 30948 BC 

The photometry of the BD companion is extrapolated as 
a relative measurement to the primary, which is described 
in Sect. |Z2] From the data obtained with the CAMELOT 
instrument, we have measured the 7-band magnitude of 
HD 130948 A to be 5.19+0.03. 

The estimation of the contrast between HD 130948 and the BD 
companion was obtain by relative comparison with a set of 
25 fake companions simulated with comparable separation and 
brightness in an image model. We compared both the peak val- 
ues of the companions and the fluxes summed in a small aper- 
ture of 3-pixels in radius. 

In order to evaluate the effect of the image post-processing on 
the BD peak value, the contrast of the fake companions in the 
image model is spanned into 50 values ranging from 5xl0" 4 to 
10~ 5 . Each of these 50 models is convolved with our observed 
PSF of HD 130948 (Image 1) and then filtered, which permits 



4. Discussion 

4.1. Prospects for the technique 

In our present work, we apply the wavelet filtering step to im- 
ages that do not benefit from AO-assisted or coronographic 
techniques. We anticipate that the assistance of, at least, a mod- 
erate AO correction able to reduce the halo contribution would 
improve the achievable contrast at shorter radii in comparison 
to what we have obtained so far. Hence, we plan to further ex- 
plore the potential of our approach with observations using the 
4.2-m WHT and its AO-system NAOMI. A fundamental aspect 
we investigated here to illustrate the potential of the technique 
is the expected sensitivity limit as a function of the distance 
to the central star. We have selected reasonably clean regions 
(i.e. unaffected by optical artefacts) of 20x20 pixels distributed 
over more than 10 azimuthal positions, centered at 3", 2", 1", 
respectively, from the center of Image 2. We have then esti- 
mated an average of the spatial standard deviation as a function 
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of the number of co-added cubes. The result is illustrated in 
Fig. 3] with the evolution of normalized SNR as a function of 
the number of co-added cubes. The bottom plot shows a reg- 
ular increase of the normalized SNR following a N~ 0A power 
law, where N is the number of co-added cubes. Closer to the 
star, the plots at 2" and 1" obtained with the current data in- 
dicate that a relatively high SNR is quickly obtained already 
after ~10 cubes, while only slowly increasing afterwards com- 
pared to the 3". Although it is difficult to derive a clear and 
decisive trend simply with the current data, these plots together 
with the filtered images presented in this work suggest that as 
close as l"-2" the sensitivity is limited by residual speckles, 
which dominate the image as for traditional AO. Hence, in- 
creasing the number of cubes under this regime does not pro- 
vide a significant gain. This effect should become limited to 
smaller separation with telescopes larger than the 2.5-m NOT. 
The small drop in the top plot of Fig.|4]around N-35 is possibly 
due to an unexpected increase of the speckle noise contribution 
in the corresponding images. Being increasingly dominated by 
speckle noise at short distance from the star, a further study of 
the speckle pattern statistics would be required, which is nev- 
ertheless not within the objectives of the current paper. 

4.2. Colour 

We place our colour measurement in the context of L- type 
brown dwarfs available photometry dLiebert & Gizis 2006). In 
Fig. [5] we reproduce a plot from these last authors giving the 
Ic-Jmko colour versus spectral type for a sample of late M 




number of cubes 

Fig. 4. Normalized SNR as a function of the number of co- 
added cubes used for the final post-processed image. The 
curves were obtained by computing the spatial rms in a 20x20 
pixels box at 3", 2" and 1" from the primary. As a reminder, 
in the present data one cube is composed of 1000 images 
with 30 ms exposure from which are selected the 30% "best" 
(brightest speckle) images. We overplotted the power-law that 
best reproduce the experimental data for 3" separation, namely 
a N~ 4 law (dotted line). The x-axis starts at value 0, and each 
sub-interval represents 5 cubes. 




Log[ peak contrast in the model image ] 

Fig. 3. Evaluation of the impact of the image post-processing 
on the peak value of a nearby companion. The solid curve 
with error bars shows, in the y-axis, the experimental mean 
and rms values of the measured peak contrast after image post- 
processed for various simulated companions contrast at a dis- 
tance comparable to HD 1 30948 BC as an input. The x-axis 
shows the input pre-determined model contrast. The dash line 
is the y-x curve and represents for comparison the contrast val- 
ues ideally unaffected by the post-processing stage. 



and L dwarfs, after conversion from the 2MASS to the MKO 
photometric system. Our first measurement of the integrated 
I c magnitude of the HD 130948 BD companion allows us to 
compare it with existing data. We assume for HD 1 30948 BC 
a spectral type L2+ 2 inferred by near-infrared spectroscopy by 



Potter et all ( 2002 ). The I- J colour me an and dispersion in the 



estimated spectral type interval for the ILiebert & Gizis (2006) 
sample equal 3.58+0.17. Within the +2 sub-class uncertainty, 
the given sample shows only one L dwarfs over 21 objects 
with bluer colour (7-7=3.19) than HD 130948 BC, and classi- 
fied L0. Hence, this suggests that our object clearly lies in the 
blue region of the colour distribution for this spectral range. 
This trend becomes even stron ger if we place our measurement 
in the context of the work by iGoto et al. (2002), who classi- 
fied HD 1 30948 BC as L4+1. A similar behavior is also found 
when placing ou r colour measurem ent i n the I-J versus spe c- 
tral type plots in iDahn et all J2002h and iBihain et all feoioh , 
where typical colours ~ 3.4-3.6 are found for L2 dwarfs. The 
field dwarf 2MASSWJ1841086+311727 (L4-type, 7 c -7«3.2, 
not plotted in Fig. [5} is another known case of spectroscopically 
reported as an L4 dwarf, although it appears too blue even for 
an early L (Kirkp atrick et alJ 2000). Ev entually, our BD objec t 
was placed in the J-K versus I-J plot oflLiebert & Giz is (2006) 
using J-K~l.6 measured by Bih ain et alJ ( I2010I) . In this plot, 
HD 130948 BC appears consistent with an L2 spectral type, al- 
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Spectral type 

Fig. 5. 7 r -7iviKn colou r versus spectral type reproduced from 



Liebert & Gizi s (2006) for late M and L dwarfs. Spectral types 
5 to 15 corresponds to M5 to L5. Error bars in the photometry 
represent the measurement dispersion (rms). Our measurement 
at the 1-cr confidence level of HD 1 30948 BC colour is shown 
as open circle. 



tional technique investigated in this paper. Alternatively, stable 
differences in the surface cloud coverage or in the atmosphere 
dust properties could alter the measured I-J index. 
An additional hypothesis, related to the presence of dust in 
the neighborhood of HD 1 30948 BC system, could be formu- 
lated as the effect of stellar radiation scattered by a residual 
of dusty disk around the BD companion, with an increased 
efficiency of the phenomenon towards optical wavelengths 
( Kal as et al.l 120071) . Using modeling tools such as the debris 
disk ra diative transfer simulator^ (DPS. IWolf & Hillenbrand 
(2005J)), we were able to assess that a small-size "circum- 
brown dwarf" disk, with a ~10~ 7 M ffi mass and composed of 
small (~0.2/im) grains confined in an ~2AU region in ra- 
dius around the BC companion, would produce a I-J varia- 
tion of ~0.2 mag towards bluer colours. However, the effec- 
tive extent and geometry of such a disk can clearly not be 
constrained with these current data. Further analytical and ob- 
servational tests (e.g. focusing on searches for small mid-IR 
excess for this object) should indicate us if such a compact 
dusty disk would survive in a relatively old multiple system, 
and how the bina ry nature of HP 13 0948 BC, with a -2.2 AU 
semi-major axis (IPupuv et al.l 120091) . or even the interaction 
with the primary would dynamical sc ulpt the dust distribu- 
tion i nto a circum-(substellar)bi narv (jArtvmowicz & Lubow 
11994 and/or circum-(sub)stellar (Mathieul l2007l) component. 
Note that, whether this hypothesis is plausible or not, the pres- 
ence of a disk would not be uncommon in sub stellar objects 
dZapatero Osorio et al.lE007HLuhman et al.ll2009l) . 



though lying slightly off the bulk of the L2 BO distribution, 
suggesting that the object may be peculiar in that diagram as 
well. Considering this trend, we intend to provide some qual- 
itative explanations to the relative blue colour of our object, 
keeping in mind that further observations would certainly help 
investigating them. 



As suggested by Mar lev et al.l (12002b . a possible reason is that 
for late L dwarfs, the shorter wavelength bands are influenced 
more by the wings of potassium and other alkali gas opacities. 
Also, the metal li city o f HD 130948 was measured by 
Valenti & Fischerl d2005l) to Fe/H=0.050+0.030, i.e. solar 



metallicity, thus presumably the same value can apply to the 
BD companion. L-type objects in the Pleiades, with metallicity 
compara ble to the Sun, do n ot present any remarkable bluer I-J 
colours jBihain et al.ll2010h . hence metallicity can probably be 
ruled out as a cause of the observed effect. 
Low surface gravity could also result in a bluer I-J colour due 
to a weaker potassium line (Kirkpatric k et al .1 120081) . Since our 
measured I-J colour for HD 130948 is slightly smaller than 
the average I-J col our in the ~120Myr Pleiades (cf. Fig. 5 in 
Bihai n et al .1 (1201 Oh l. we may have an indication for an even 
younger age of HD 130948 in compa r ison t o the 0.1-1 Gyr 
range of age proposed by IPupuv et al . (2009). Also, weather- 
like phenomena in the BD a tmosphere dBailer- Jones & Mundt 
19991: Gold man et al.ll2008l) could cause flux time-variability 
resulting in a larger dispersion of the I-J colour. Such a case 
may be explored with repeated monitoring observations of the 
object, which could be possibly conducted using the observa- 



5. Conclusions 

This work addresses the question of high contrast in optical 
speckle imaging in the context of substellar objects, in a spec- 
tral domain where conventional AO systems present modest 
performances for long-exposure imaging. We show that the 
Lucky-imaging approach has a significant potential for the de- 
tectability of faint and close companions to bright stars, which 
can be improved with an additional stage of post-processing 
based on the image wavelet transform. We successfully applied 
this approach for the first time to the system HD 130948 with 
a clear detection in the I c filter of the BD binary companion 
HD 130948 BC at -2.5" using the 2.5-m NOT telescope, and 
in the absence of any coronographic system. The BD compan- 
ion is spatially unresolved because of a too small physical sepa- 
ration (30 mas) at the time of our observations, which is below 
the telescope resolution at this wavelength. The relative con- 
trast at this distance is estimated to Al = 11.30 (~3.0xlO -5 ). At 
1", the point-source detectability is estimated at A/~10 in the 
neighborhood of a bright star like HD 130948 A, making the 
proposed approach a powerful technique to obtain high resolu- 
tion and high-contrast photometry even with small class tele- 
scopes. Considering these results, we plan hopefully a positive 
extension of this a pproach to larg er and AO-assisted medium 



class telescopes (lLaw et al.l |2009) or eventually to corono 



graphic imaging system. Alternatives to substellar programs 
could also be considered for science cases where optical imag- 
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ing is an asset, as for instance in the case of the search for white 
dwarf companions to bright stars. 

Appendix A: Comparison and details on different 
image processing algorithms 

In this work, three different algorithms to enhance the detec- 
tion of HD 1 30948 BC. For all of the three options, the filters 
are applied recursively in order to improve the companion de- 
tection (~3 iterations, see below). These algorithms are: 

★ a median box filter, 7=/-median(/,N), where N is the size of 
the median box filter, followed by a convolution with Gaussian 
kernel to filter out pixel noise. 

★ an unsharp mask filter, I=I-I*g, where g is a Gaussian ker- 
nel, followed by pixel noise filtering. 

★ a wavelet mask filter resulting from the wavelet transform 
of the image. As a multi-resolution decomposition of the origi- 
nal image, this transform provides a set of filtered images from 
which can be selected the best trade-off between the spatial res- 
olution and the SNR. The basic algorithm used in this work is 
largely inspired from the wavelet "a trous" procedure written 
in C/Yorick language and part of the Yeti package (see www- 
obs.univ-lyonl.fr/labo/perso/eric.thiebaut/yeti.html). We apply 
to the dataset a the function WT(a,ORDER), where WT is the 
wavelet transform described hereafter in the one-dimensional 
case and extendable to the two-dimensional case of images 
for our work. The parameter ORDER gives the maximum de- 
gree of decomposition in the multi-resolution analysis. In the 
following, i is the intermediate decomposition degree (also 
called order) and 1^/^ORDER+l. By applying the operation 
/?=WT(a,ORDER), b becomes a cube of images of the same 
size as a and with its third dimension equal to ORDER+1, 
as explained hereafter. The image b\ of intermediate order i is 
given by: 

a\ — a (A.l) 
bi — a\ - Oi+i (A.2) 
ai+i = T(Oi, SCALE = 2 i_1 ) (A.3) 
T is a convolution operation by: 
a i+ i(p)= ^0 + W + 1) x ai(p + j x SCALE) (A.4) 

where SCALE=2'~' . Integer p is the pixel number in the image, 
"SCALE" and W are, respectively, the scaling factor and the 
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Fig. A.l. Effect of the recursivity in the image filtering process. 
This serie shows our resulting filtered image of the HD 130948 
system after one (left), two (center) and three (right) iterations 
of the wavelet algorithm used in this case. 



Fig. A.2. Comparison of the different algorithms used in this 
study. The intensity scale has been optimized to detect the 
fake companions. 10 individual companions have been simu- 
lated for the tests. Only 5 are displayed here. All the presented 
images have undergone 3 iterations of the tested filtering al- 
gorithms. Left: wavelet algorithm used in this work. Upper- 
middle: 20x20 pixels Gaussian kernel filtering followed by 
3x3 Gaussian smoothing to remove pixel noise. Lower-middle: 
5x5 pixels Gaussian kernel filtering followed by 3x3 Gaussian 
smoothing to remove pixel noise. Upper-right: 20x20 pixels 
median filtering followed by 3x3 Gaussian smoothing to re- 
move pixel noise. Lower-right: 5x5 pixels median filtering fol- 
lowed by 3x3 Gaussian smoothing to remove pixel noise. 



half-size of the wavelet kernel k (i.e. k has a length of 2xW+l 
pixels in the one-dimensional case). From cube b, we selected 
one filtered image, at the optimal order which represents the 
best trade-off between the necessary spatial resolution to de- 
tect the companion and a good SNR. If required, the selected 
filtered image becomes the input image a and the filtering pro- 
cess is run recursively, keeping the same value of the ORDER 
parameter and eventually selecting the same optimal interme- 
diate order. The reason why this process is run recursively - 
or iteratively - is to achieve a higher level of suppression of 
remaining "low spatial frequency structures" in the image, as 
illustrated in the sequence of Fig. lATl We verified empirically 
that there is generally no significant gain after a maximum of 3 
iterations. 

The wavelet kernel is support-limited and has been em- 
pirically optimized for this work to a 5x5 pixels square 
(i.e. W-2) with a two-dimensional distribution given by 
[[1,1, 1,1,1], [1,3,3,3,1], [1,3,8,3,1], [1,3,3,3,1], [1,1, 1,1,1]]. 

In order to evaluate the effects of each algorithm, we 
have created 10 fake companions with similar properties to 
HD 1 30948 BC and radially distributed around HD 130948 A. 
As an additional test, this permitted us to discard possible al- 
gorithm artefacts that would have led to a false detection of 
HD1 30948 BC. 

The imaging results are shown in Fig. IA.2l We initially imple- 
mented the two first algorithms using a kernel with a FWHM of 
4-5 times the original PSF FWHM. In both cases, this has led 
to an excessive smoothing of the image preventing us from de- 
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tecting all of the companions and in particular HD 1 30948 BC 
(see Fig. IA.21 upper-middle and upper-right panel). Reducing 
the convolution kernel to a size comparable with the PSF 
FWHM improved significantly the result, although the median 
box filter does not reach a satisfactory level. On the contrary, 
the iterative unsharp mask filtering using a smaller kernel per- 
mits us to detect all the companions, included HD 130948 BC. 
This suggests that, apart from the multi-resolution decomposi- 
tion - which is useful to identify a close-to optimal kernel size - 
unsharp mask filtering is almost identical in terms of point-like 
source detectability to the wavelet filtering algorithm described 
above. 
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